Visibility of graphene flakes on a dielectric substrate 
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We model the optical visibility of monolayer and bilayer graphene deposited on a Si02/Si sub- 
strate or thermally annealed on the surface of SiC. Visibility is much stonger in reflection than in 
transmission, reaching the optimum conditions when the bare substrate transmits light resonantly. 
In the optical range of frequencies a bilayer is approximately twice as visible as a monolayer thereby 
making the two types of graphene distinguishable from each other. 



Monolayer graphene is a single two-dimensional hon- 
eycomb lattice of carbon atoms. Although the first 
graphene-based structures were only recently fabricated 
[TJ they have quickly become the subject of an exten- 
sive research effort [2] |3j 3]. Monolayer graphene is a 
zero-gap semiconductor with a Dirac-like dispersion of 
chiral quasiparticles near the K points of the hexago- 
nal first Brillouin zone [5.. Bilayer graphene is a pair of 
graphene sheets with the Bernal (AB) stacking arrange- 
ment. In the low-energy spectrum of this material [6 
the conduction and valence bands both consist of two 
quadratic branches split by the inter-layer coupling 71. 
Measurements of the quantum Hall effect [TJ |2] [7] and 
ARPES experiments [8] have confirmed that these are 
the low-energy band structures of these materials. 

The widespread microcleavage technique used to fabri- 
cate graphene-based devices requires a visual inspection 
of the substrate [TJ to find flakes of one or two layers 
thickness. In this Letter, we aim to determine the opti- 
mum conditions for making these flakes optically visible 
when they are deposited on various substrates. The pa- 
rameters at one's disposal (see Fig. [TJ are the frequency 
u>, angle a. and aperture 5a of the focused incident radi- 
ation, as well as the thicknesses of the various layers of 
the underlying dielectric materials. 

Below we calculate the reflection of non-polarized in- 
cident light taking the geometry of the substrate into 
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Figure 1: Geometrical configuration for detection of graphene 
on a substrate. A light beam is focused on a small spot which 
is scanned along the surface. The calculations presented be- 
low show how to optimise conditions for visibility of atomi- 
cally thin graphitic flakes. 



account with suitable boundary conditions at each of 
the interfaces between materials, appropriate frequency- 
dependent dielectric functions e(u>) for each layer, and 
fj, = 1, Throughout the calculation, we use the data [9] 
available in the existing literature for the dispersion of 
the permittivity of silicon [10], silicon oxide [11] and sili- 
con carbide [12]. With reference to Fig. [TJ we analyze the 
reflection R of light from a substrate with a flake on it 
and compare this to the reflection R of a bare (graphene- 
free) substrate. The optical visibility of a flake is then 
determined as the contrast between two such parts of the 
sample studied using a monochromatic light source: 

V R = (R-R )/R . (1) 

The scattering of light is analyzed using the electro- 
magnetic wave equations in vacuum and dielectric media 
and the standard boundary conditions at interfaces be- 
tween different materials, 

E\= e\, E>i =D^, B,\-Bl= (t(uj)E x n. (2) 

The superscripts || and _L stand for the components of the 
field parallel and perpendicular to the interface respec- 
tively, n is the unit vector normal to the interface, the 
subscript 1 (2) denotes the field below (above) the inter- 
face, and <j(oj) is the frequency-dependent conductivity 
of a graphene flake and D — e(u>)E. One more bound- 
ary condition (on the perpendicular components of H) 
duplicates Snell's law. 

Having in mind an optical setup used to locate a small 
flake, we consider a beam of light focused by a lens, so 
that the light in the beam arrives at the substrate sur- 
face with some aperture 5a (see Fig. [T}. Therefore the 
measurable reflectance to be used in Eq. (JTJ is 

R(a,5a) = J (Kl n R(k)P(k), (3) 

where P(k) characterises the spread of the beam over 
the solid angle of the aperture 5a around a, k — 
^(sin a, 0, — cos a) is the wave vector of the incident ray 

of light, and R(k) is the reflection coefficient for a plane 
wave with this wave vector. Below we assume that the 
beam is equally dense at all angles within an aperture of 
5a around a. 

To describe the conductivity of graphene, we follow the 
method used in Refs |13l IT3] taking into account the split 
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bands formed in the bilayer [6]. At low temperatures the 
result for the monolayer which takes into account the 
transition between the valence and conduction bands in 
the Dirac spectrum is <j\ = e 2 /Ah (with a negligible imag- 
inary part) [13]. This corresponds [15] to the absorption 
coefficient g = Ana/c which gives g\ — ire 2 /he » 2.5%. 
For the bilayer, there are four possible inter-band transi- 
tions, reflected by its conductivity, 



0(0 - 2) 




(4) 



Here ft — Huj/ji is the frequency written in units of the 
inter-layer coupling and 8(x) — [1 + sgn(a;)]/2. The real 
part of this function has a discontinuity at huj — 71 ss 
0.4eV and a cusp at hui — 271. These correspond to the 
activation (at zero temperature) of the interband transi- 
tions between low-energy bands and split band, and the 
two split bands respectively. The imaginary part of 02 
shows a divergency at hui = 71, leading to an enhanced 
reflectance of the bilayer at this frequency. 

For non-polarized light arriving at the incidence an- 
gle a to the sample depicted on the right-hand side of 
Fig. [TJwith graphene deposited on the top surface, the 
reflectance is 
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In this result the first term represents reflection of radi- 
ation polarized so that the electric field is perpendicular 
to the plane of incidence, the second term to radiation 
polarised so that the electric field is parallel to the plane 
of incidence, and 

A = — \fs~d cos ad cos X<i + iy/sbcos a^sin Xd, 
B = i\fsd cos a d cos Xd — \feb cos a^ cos X^, 
C = — i^/iEd cos a dB sin X s + y^cos a s Acos X s , 
D = \fe~d cos ctdB cos X s — zy^cos a s AsinX s ; 

A' = y/Ii, cos ad cos Xd — iy/sa cos a b sin Xd, 
B = y/ed cos ab cos Xd — iy/ei cos ad sin Xd, 
C = yfed cos a s A' cos X s ~ i^/Fs cos adB' sin X s , 
D' = — i y/id cos a s A 1 sin X s + -y/ejeos adB' cos X s . 



Here X s = y^fcs cos a s , Xd — y/edkd cos ad, sina^ = 
sina/y 7 ?^, sina s = sina/^/e^ and sincy = sin a / y/Fd- 
The and a is determined by the direction of the wave 



vector of the incident plane wave, see Fig. [T] To model a 
finite slab of silicon of width d with a silicon oxide layer 
of width s on top, we substitute Ed — £^i(u)), e s — £si0 2 i 
£5 = 1, and the quantity Rq is found by replacing a = 
in these expressions. To evaluate the visibility Vr, the 
integral in Eq. ((3I must be taken for R and R$ using Eq. 
d5l 
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Figure 2: Visibility in (a) reflectance Vr and (b) transmit- 
tance Vt for graphene on a silicon wafer with a 300nm oxide 
layer for varying substrate thickness and frequency of radia- 
tion. Note the difference in the scales of (a) and (b). Here we 
take an aperture of 8a — 10° and a = 20°. 

Figure [2] illustrates the visibility of mono- and bilayer 
flakes on a Si substrate of widths 0.5/im < d < 1.5/im and 
a 300nm Si0 2 layer (see Fig. [TJ for light with 0.3eV < 
hui < 2.5eV arriving with aperture da = 10° around a 
20 * . The rapid oscillations of the visibility in this plot are 
caused by the resonant condition of the Si layer. When 
this layer is strongly transmitting (that is, when cosXd ~ 
0), the visibility is at its highest. This fine structure is 
modulated by the corresponding resonance condition in 
the oxide which is responsible for the 'bands' which lie 
across the plots in Fig. [2] The condition for maximum 
transmission through the oxide is cosX s w which leads 
to 



uj 1=3 c(n + |)7t/ (sy/e^cos a s ) 



(6) 



where n is an integer. The wave vector of the light in 
the slab is of the order of an inverse micron, so the res- 
onant conditions are closely spaced on the length scale 
of the substrate thickness. The visibility of a bilayer 
flake is higher than the visibility of a monolayer for 
Hui > 7i w 0.4eV because the conductivity of the bi- 
layer is essentially twice as large as the conductivity of 
the monolayer in this energy range. Additionally, the 
divergency in the imaginary part of the bilayer conduc- 
tivity at huj — 71 rj 0.4eV causes a stronger reflection 
and hence a larger visibility. Also we have calculated the 
transmittance T of the sample, and the corresponding 
visiblity V T = (T - T )/T is shown in Fig. gb) where 
the same resonant structure appears, but is at least ten 
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Figure 3: Frequency dependence of visibility Vr of graphene 
on (a) an infinite silicon slab with a thin oxide layer of width 
300nm, and (b) a silicon carbide slab of width 1/im. In both 
plots we use aperture 8a = 10°. 



times weaker than the visibility in reflectance. 

We find that the visibility of graphene in reflectance 
is futher enhanced by using a thick (semi-infinite) sub- 
strate with a sizeable oxide layer on its surface, in agree- 
ment with a recent experimental observation [16 . Fig- 
ure |3ja) shows the visibility of graphene deposited on a 
semi-infinite slab of silicon j9] with a 300nm SiC>2 layer. 
In this case the analytical expression for the reflectance 
of a plane wave with wave vector k = — (sin a, 0, — cos a) 
can be found by substituting £& = esi(of), e s — £sio 2 



and d = into Eq Q5B . As before, both R(k) and Ro(k) 
(which is determined from this equation with er = 0), 
must be substituted in Eq. ([3} before the visibility is 
evaluated. In the plots in Fig. |3^a), the main features 
are the very strong reflectance of the graphene flake at 
huj w 0.5eV and hw w 1.6eV. These are due to the stand- 
ing wave resonances in the oxide layer at the condition in 
Eq. (|6}. In Fig. |3ja) the peak in visibility at Tiw « 0.5eV, 
(n = 0) corresponds to the first resonance in the oxide 
layer and the peak at fuv w 1.6eV (n = 1) to the second 
resonance. The factor of 2 difference between the bilayer 
and monolayer conductivities at hut ^> 71 and the diver- 
gence in the imaginary part of 02 (^0 at Huj = 71 0.4eV 
are manifested in the visibility. 

Besides being produced using the microcleavage tech- 
nique, ultra-thin graphitic films can also be grown by 
thermal annealing of SiC wafers [HE]- The reflectance 
for this configuration can be found by substituting [9 
d = 0, Eb = 1 and e s = esic m Eq. Plots of the 

visibility defined by this function are shown in Fig. [3](b) . 
The standing wave resonance in the substrate is again 
the main factor for the visibility of graphene, though it 
is weaker for a SiC slab than for the SiC^/Si substrates. 

In conclusion, we have found that graphene is much 
more visible in reflection than in transmission and that 
the resonance condition of the substrate is the dominat- 
ing factor in determining its visibility. For optimum visi- 
bility the wavelength of monochromatic light used should 
be selected using Eq. (16} , and for the visible frequency 
range where {a% ~ 2a\) a bilayer is clearly distinguish- 
able from a monolayer. 
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